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^ : Abstract 

' We continue our study of patch effect (PE) for two close cylindrical conductors with parallel 

axes, slightly shifted against each other in the radial and by any length in the axial direction, 

o'! r 

Q" 1 started in where the potential and energy in the gap were calculated to the second order in 

the small transverse shift, and to lowest order in the gap to cylinder radius ratio. Based on these 
results, here we derive and analyze PE force. It consists of three parts: the usual capacitor force 
due to the uniform potential difference, the one from the interaction between the voltage patches 



and the uniform voltage difference, and the force due to patch interaction, entirely independent of 
the uniform voltage. General formulas for these forces are found, and their general properties are 
described. A convenient model of a localized patch is then suggested that allows us to calculate all 
the forces in a closed elementary form. Using this, a detailed analysis of the patch interaction for 
one pair of patches is carried out, and the dependence of forces on the patch parameters (width and 
strength) and their mutual position is examined. We also give various estimates of the axial patch 
effect force important for the Satellite Test of the Equivalence Principle (STEP), and recommend 
intensive pre-flight simulations employing our results. 
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I. INTRODUCTION 

Electrostatic patch effect |2j is a nonuniform potential distribution on the surface of a 
metal first examined theoretically in paper [3J, where the PE force in a plane capacitor with 
varying boundary voltages has been calculated. This analysis wasparticularly motivated by 
the LISA space experiment to detect gravitational waves (see c.f. [4J). PE can similarly affect 
the accuracy of any other precision measurement if its set-up includes conducting surfaces in 
a closed proximity to each other. PE torques turned out one of the two major difficulties (s| 
Q] in the analysis of data from Gravity Probe B (GP-B) Relativity Science Mission (2004- 
2005 flight), that measured relativistic drift of a gyroscope predicted by Einstein's general 
relativity This required theoretical calculation of PE torques jsj for the case of two 
concentric spherical conductors. 

Here we continue studying PE in cylindrical geometry that we started in paper hence- 
forth referred to as CPEI, for "'Cylindrical Patch Effect"'. It is strongly motivated by the 
experimental configuration of STEP 10h13|. where each test mass and its superconducting 
magnetic bearing is a pair of approximately coaxial conducting cylinders (see section IVI Al 
for more details). The goal of STEP is the precise (1 part in 10 18 ) measurement of the 
relative axial acceleration of a pair of coaxial test masses, so the importance of properly ac- 
counting for PE forces is evident. Notably, the axial and transverse force we find is inversely 
proportional to the gap and its square, respectively, exactly as in the plane capacitor j^|. 

We determine the PE forces between the two slightly shifted cylinders with parallel axes 
by the energy conservation argument: a small shift, f °, of one of the conductors relative to 
the other causes an electrostatic force given by (see, for instance, [yj]) 



(i) 

f o = o 



where W(r°) is the electrostatic energy as a function of the shift. The latter was found in 
CPEI to the second order in the small transverse shift, po *C d, where d is the gap between 
the cylinders in the coaxial position, so the force to the first order in pn is found below. For 
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121 ]. the gap is much 



typical experimental conditions, such as the STEP configuration 
smaller than either of the cylinders' radii; thus two small parameters are actually involved in 
the problem, p /d <C 1, d/a ~ d/b <C 1. While justifying the model of infinite cylinders, 
this also allows for a significant simplification of the results to 1. o. in d/a. 
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In the next section we summarize the results from CPEI needed for the force calculation. 
Based on this, we calculate PE forces in section IHIt and illucidate their general properties. 
In section IIVI we introduce a convenient model of a localized patch potential allowing one 
to find simple closed-form expressions for the forces. Section |V] contains a detailed analysis 
of PE forces when a single localized patch described by our model is sitting at each of 
the cylinder boundaries. In section I VI I estimates of the axial patch force for the STEP 
experiment set-up are given. The details of calculations, in places rather complicated and 
cumbersome, are found in the appendix. 



II. SUMMARY OF RESULTS FROM CPEI 



We use Cartesian and cylindrical coordinates in two frames of the inner and outer cylin- 
ders as shown in fig. [TJ In the inner, or 'primed', frame the position of a point is given 
by the vector radius r , and Cartesian coordinates {x ,y',z'} or cylindrical coordinates 
{//, ip , z }. In the outer, or 'unprimed', frame the corresponding quantities are r, {x, z}, 
{p, ip, z}. Frame origins are separated by r*°, hence the coordinates are related by 

f =f+f°; x =x + x°, y =y + y°, z = z + z° . (2) 

We also use an alternative writing x = x®, y = x 2 °, x = xi, y' = x 2 , etc. 

The surfaces of the inner and outer cylinders are p' = a and p = b, respectively, with 
d = b — a. They carry arbitrary voltage distributions given by the conditions 



= V-+H(<p,z), (3) 

P =b 



= G(<p,z), $ 

p =a 

where $ is the electrostatic potential, and V~ is the uniform potential difference: all the 
voltages in the problem are counted from the uniform voltage of the inner cylinder taken 
as zero. The non-uniform potentials (patch patterns) are described by arbitrary smooth 
enough functions G(ip' , z ) and H(<p, z), whose local nature is emphasized by requiring 

2ir oo 2ir oo 

||G|| 2 = y J dtp'dz' \G{^',z')\ 2 < oo, \\H\\ 2 = J J d^dz\H{ v ,z)\ 2 < oo . (4) 

— oo — oo 

For any function u(<p, z) satisfying the square integrability condition we denote its Fourier 
coefficient u n (k): 

2lT oo 

1 /"°° _ i r r 

u(<p,z) = — / dk £ u n (k)e^ kz+ ^\ u n (k) = — / / d<pdzu(cp,z)e-« kz+n ^; (5) 



in particular, Fourier coefficients of the patch voltages G(ip ', z ') and H(ip, z) are G n (k) and 
H n (k), respectively. Since these functions are real, their Fourier coefficients satisfy 

G n (k) = G*_ n (-k), H n {k) = H*_ n (-k) ; (6) 

here and elsewhere the star denotes complex conjugation. For any two squarely integrable 
functions u(ip,z) and v((p,z) the useful Parceval identity holds: 

27T CO 



/ / d(pdzu((f, z)v*(ip, z) = I dk ^2 u n (k)v*(k) . (7) 



-oo 



In the case u = v the identity ([7j) shows that the squared norm, \\u\\ 2 , of a function u is 
equal to the squared norm of its Fourier coefficient u n (k). 

As shown in CPEI, section III, electrostatic energy in the gap between the two cylinders, 
as a function of their mutual shift r consists of three parts, 

W(f°) = W u {f°) + W int (f °) + W p {f °) , (8) 

where the first one is due to the uniform potential difference, the second results from the 
interaction between the uniform and patch voltages, and the third one is the energy of 
the patch interaction. Each of these contributions was found in CPEI in the form of an 
expansion in the small transverse shift po to quadratic order in the small po/d ratio, with 
the coefficients depending generally on the axial shift zq\ 

W A (f°) = W A (z°) + Wf(z°) (xl/d) + W*(z°) (xl/d) (x°Jd) + O [(po/rf) 3 ] ; 

A — u, int, p . (9) 

The explicit coefficients for each kind of energy follow immediately. 

The uniform potential energy is, naturally, proportional to the length of the capacitor, 
same as its expansion coefficients which, for \z\ < L, are 

W U (L) = 2-nLt,- d (V~)\ W;{L) = 0, W^(L) = nLe ^ (vf 5, u . (10) 

Here and everywhere else we adopt the summation rule over repeated Greek indeces p, u, 
etc.: the summation over them runs from 1 to 2, corresponding to transverse coordinates in 
the cylinder cross-section. The first order term vanishes as it should be due to symmetry 
(otherwise there would be a non-zero transverse force in a perfectly symmetric configuration 
of coaxial conductors under uniform potentials). 



The expansion coefficients for the interaction energy are: 



W, 



-2ne -V- [(G (0) - H (0))} , 
+A7re ^V-^\cUG 1 (0)-H 1 (0)) 



int 

" "^d >' 



(11) 



Wj* = -^- d V ~ [ c >t (G a (0) - H 2 (0))} + (V/4) (G o (0) - H (0))} 



where cf = 0.5, c 2 = 0.52, 3?(-) denotes the real part of (•); recall also that G n (k) and H n (k) 
are the Fourier coefficients of the patch voltages (jBJ). Same as the uniform part of energy, 
the interaction one does not depend on the axial shift, because the longitudinal shifting of 
an electorde with the uniform potential does not change the actual charge configuration. 
Finally, the patch energy coefficients are found to be: 



w p 

ft 



2d 



/oo 00 n 
dk £ \G n (k)e lkz -H n (k)f 
-00 _„ 



e a 



, 1 dk ^ \ c t (Gn(k) e ~ lkz ° ~ H*{k) \ (G n+1 (k)e tkzU - H n+ i(k) 

CL J — oo - 



(12) 



n=— 00 
oo 



W*, = / dk Y: {S, v /2\G n (k)e^° - H n (k)\ 2 + 
za J-00 



n=— oo 



-tkz° 



K(k)) (G n+2 (ky kz ° - H n+2 (k))}} . 



These quantities do depend on the axial shift, z°. Combining expressions flS} and (J2J) we can 
also write for the total energy: 



W^ ) = W^ 1 + Wf ' + (A £ = 0, /x, /ii/ . 
Note that all parts of energy are given to 1. o. in d/a. 



(13) 
(14) 



III. PATCH EFFECT FORCES 



As explained in the introduction, the force is found by the formulas ([T]) and ([TBI : 

dW(f°) 



dxl 



W^) + 2W llu (z») {xl/d) + O {{p /d) 



^(0) + W'(0) z° + 2W^(0) {x Q Jd) + O ((r /d) 



1,2 



(15) 



\f °| = 



vV) 2 + (2/ ) 2 + (*°) S 
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for the transverse force components, and 

dW(f°) 



dz° 



- W' {z°) + W;(z°) (xl/d) + O {{po/df) 
- <(0) + <(0) z° + <(0) (xl/d) + O ((r /rf) 2 )] , (16) 

for the axial force (here and everywhere else primes denote the derivatives in z°). Same as 
energy, the force consists, of course of three parts, which we study below one by one. 



A. The Force due to the Uniform Potential Difference 

Using the formula (fT5l) for the transverse force and expansion (Q with the coefficients 
(HUD we obtain (as usual, F? = F% = F^): 

F: = -2nLe -^(V-) 2 (x°/d), F; = -27rLe ^ (V') 2 (y°/d) , F z w = . (17) 

Zero axial force, obvious by symmetry, is formally due to the independence of the energy on 
the axial shift, z°. Recall that 2L is the cylinder length, so the force per unit length of the 
cylinders that is finite. 



B. The Patch and Uniform Potential Interaction Force 

To get it, we combine formulas f lT5|) . fl9]) and f[TT|) ; after some simplifying transformations 
the result becomes [9(-) is the imaginary part of (•)]: 

F™ = ~^o^V- {3ft [G a (0) - 11,(0)} - 

3ft [(G 2 (0) - H 2 (0)) + (G o (0) - H (0))] (x°/d) + 3 [G a (0) - H 2 (0)] (y°/d) } ; (18) 
= -2ne -^V- {-9 [Gi(0) - #i(0)] + 9 [G 2 (0) - H 2 (0)} (x°/d) + 
3ft [(G 2 (0) - H 2 (0)) - (G (0) - # (0))] (y°/d) } ; 

F l z nt = . 

The axial force vanishes again, by symmetry. 



C. Forces due to the Patch Interaction 

By the formula (j!5p and the energy expansion (JHJ) with the coefficients (1121) we find: 

la j-oo n= _ 00 
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G* n {k)e 



-tkz° 



H* n (k)) [G n+1 (ky kz °- H n+1 (k) 



G n (ky kz ° - H n (k)\ 2 + 3ft [(G;(fc)e- 4fc2 ° - H*(k)) (G n+2 (k)e ikz ° - # n+2 (*0)]] (*7<*) + 

3 [(G;(A;)e- fez0 - H* n (k)) (G n+2 {k)e lkz ° - H n+2 (k))} (y°/d)} ; (19) 

en r oo 

^ P = - / dfe £ {3 [te(^e"^° - (G n+1 {k)e lkz ° - H n+1 (k)\ 

la j-oo n= _oo 

G* n (k)e- tkz ° - #*(*;)) (G n+2 (A;)e^ - H n+2 (k))] (x°/d) + 
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G n (k)e tkz ° - H n (k)\ 2 - 3ft [(G;(A;)e-^ - H*(kfj (G n+2 (k)e tkz ° - tf n+2 (A;) 
To get the axial force, we use formulas ([9]), (|T2|) . and (|T6l) : 



V")} • 



+ 
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F * = - -T- / E ( 9 [fcGn(fc)#n(fcK 

« n=-oo 

G* n (k)H n+1 (k)e-* kz ° - H* n (k)G n+1 (k)e lkz °} k (x°/d) + 
~G:{k)H n+1 (k)e~ lkz ° - H* n (k)G n+1 (ky kz °] k (y°/d)\ . 



(20) 
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Formulas ffl9|) and (120]) provide the force due to patches to linear order in a small trans- 
verse shift for an arbitrary axial displacement. In many cases, such as the STEP set-up 
described in sec. IVI| the axial shift is also small, and PE forces to linear order in all the 
shifts are only needed. We obtain these expressions by replacing exp (±tkz°) with the two 
terms of its Maclaurin expansion, assuming that all the arising integrals in k converge: 



x 2d 2 



/OO 00 
dk £ {9ft [{G* n {k) - H*(k)) (G n+1 (k) - H n+1 (k))] 
-oo __ „ 



\G n (k) - H n (k)\ 2 + 3ft [(G* n (k) - H*(k)) (G n+2 (k) - H n+2 (k))}] (x°/d) + 
9 [(G* n (k) - H* n (k)) (G n+2 (k) - H n+2 (k))\ (y°/d) + 
%(G n (k)H* n+1 (k) - H n (k)G* n+1 (k)) (kz )} ; 

F V = "IS r ^ £ & ^ G nW ~ H n( k )) ( G n + l(k) ~ H n+l (k))] - 

3 [(G* n (k) - H*(k)) (G n+2 (k) - H n+2 (k))) (x°/d) + 
\G n (k) - H n (k)\ 2 - 3ft [(G* n {k) - H* n {k)) {G n+2 (k) - H n+2 (k)))] (y°/d) + 

^(G n (k)H* n+1 (k) - H n (k)G* n+1 (k)) (kz )} ; 

F* = / dk £ [kG n (k)H:(k)} + 

« •'-oo n=-oo 

- [G* n {k)H n+1 {k) - H* n (k)G n+1 (k)} k (x°/d) + 



(21) 



3ft 



[G* n (k)H n+1 (k) - H* n (k)G n+1 (k)} k (y°/d) + U{G n {k)H* n {k))k (kz° 



Unlike the transverse shifts x° and y°, the axial shift enters here not in the ratio to the gap, 
d, but in the product with k, which is the inverse characteristic length in the axial direction. 

Note that all the forces are derived as acting on the inner cylinder, the forces on the outer 
one have the opposite sign. 



D. General Properties of Electrostatic Patch Effect Forces 

The above results allow for some general conclusions regarding the patch interaction. 

1. To lowest order, the axial patch force is inversely proportional to the gap width, the 
transverse force components go as its inverse square. 

2. Forces and shifts are directionally coupled: a transverse shift causes generally some 
axial force, and vice versa, a transverse force appears due to an axial shift. 

3. The axial force vanishes when the patches are present on one of the cylinders only, 
i.e., when either G n {k) = 0, or H n (k) = 0. 

4. The transverse force does not vanish when the patches are on one of the cylinders only. 
Moreover, expressions (|T9l) for its components can be given in terms of the non-zero patch 
potential (and not its Fourier coefficient!) using the Parceval identity (j7]). For example, in 
the case when there are no patches on the inner cylinder [V a {(f ,z ) = 0, G n {k) = 0] the 
transverse components are: 

cos (f — (1 + cos 2(f) (x° I dj — sin 2(f (y / d^j ; 
sin f — (1 — cos 2f) (y / dj — sin 2(f (x° /dj . 



n = ir 2 \ dz dfV b \f,z) 



F v = £-2 dz d W b 2 {^z) 



The transverse interaction force also does not vanish in such a case, by the formulas ( TL8l) . 

5. Uniformly charged cylinders give rise to a mutual restoring (elastic) force, so that the 
coaxial condensor is at a neutrally stable equilibrium. 

6. The interaction between patch and uniform potentials involves only the first and 
second polar angle harmonics of the patch distribution if the force is taken, as above, to 
linear order in the transverse shift. 

Some of the above conclusions might be rather obviuos or intuitively clear, however, all 
of them are now accurately established by our analysis. 
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IV. THE PATCH MODEL 



To understand better patch interaction, it is natural to examine the case when just a 
couple of patches are present, and the PE forces are described by as simple expressions as 
possible. One thus needs some convenient model of a patch as a localized deviation from the 
uniform potential described by some particular functions with just few parameters involved. 
One of them should control the patch potential, two more have to govern the spot width in 
the axial and azimuthal directions, and two more parameters specify the patch position. 

Such a patch model is desired for another reason as well. In the experiments there is 
usually no way to directly measure the patch distribution at the electrode surfaces. In- 
stead, one should infer it from some other signals, like the patch forces. However, our force 
formulas are not fit for immediate modeling, since the unknowns in them are the Fourier 
coefficients G n (k) and H n (k), with no means to estimate these functions unless properly 
parameterized. The existing experience of such parameterizations, along with the common 
sense, demonstrate clearly that only the models based on the underlying physics, rather 
than ad hoc ones, turn out efficient and work successfully. 

So, the goal of an effective patch model is to find such functions that: a) both Fourier 
coefficients G n (k) and H n (k) are found in a closed form, and b) all the series and integrals in 
the formulas (fT9]) and (120]) for the forces are computed analytically in a closed form. From 
this standpoint, separation of variables is the simplest representation: 

V(y? -f^z- z*) = V, f(z - z*) u(<p - p*) . (22) 

Here the normalizing constant V* has the dimension of a potential, and the dimensionless 
functions f(z) and u(<p) are chosen so that \f(z)\ < 1, \u((p)\ < 1; /(0) = 1, u(0) = 1. 
The center of the patch is at ip = ip#, z = z*, where the potential achieves the maximum 
magnitude V m (positive or negative). The Fourier coefficient of the function ( 12 2 p is: 

V n (k) = V, j(k)e- lkz * u n e~ m ^ = K f(k) u n e~^ +n ^ ; (23) 

1 roo j /"27T 

/(*) = _=/ dzf(z)e- ikz , u n = ^= d V u{ V )e- in *. 

V27T J-oo V27T Jo 



Successful implementation of the above requirements a) and b) is in the choice of functions 
f(z) and u(ip). The first of them is natural to choose as the Gaussian exponent in z, 



f(z) = exp 



f(k) = Azexp 



'kAzX 



(24) 
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with the Fourier coefficient a Gaussian exponent, too, and the parameter Az giving the axial 
half-width of the spot. Since the product of two Gaussians entering the force expressions 
(1T9"]) and (|2"U|) is again a Gaussian, the integrals in k there combinations of the elementary 
functions, as desired. A good choice of the second function, u(p), which is 27r-periodic, 
turns out much more difficult. Nevertheless, eventually one can come up with the following: 



u((p) = u(<p,X) 



A) 5 



1 + cos p 



(25) 



2 l-2Acos^ + A 2 ' 
where — 1 < A < 1 is controlling the width of the peak at p = (see fig. [2] below). 
Indeed, when A = —1, u(ip, —1) = 1, for A = 0, u(p,0) = 0.5(1 + cos^), and finally, when 
A — > 1 — 0, the function demonstrates a 'bounded delta-like' behavior by going to zero 
everywhere except p = 0, where the limit is unity (zero width peak). We introduce the 
azimuthal patch half-width, Aip, in an accurate way as the abscissa at the point where u 
coincides with its mean value over the whole interval, u(Ap) = u av , which gives 

cos Ay? = A, Ap = arccos A . (26) 

In a complete agreement with the above, when A — > 1 — 0, the width shrinks according as 



Ap m y 2 (1 — A) — > 0; in the opposite case A = — 1 one naturally has 2A<p = 2n. 

What makes the choice (125]) really invaluable for calculations is its Fourier coefficients 
gotten by simply expanding the function in powers of Aexp (— ip): 



u n (X) 



'2ti 



1-A 2 



A W , n^0; 



u = w (A) 



27T 



A 



(27) 



which are essentially just exponents of |n|, as in the geometric progression. Apparently, 
expressions (125]) and f l27|) satisfy our requirements a) and b), perhaps even in the simplest 
possible way. The profiles of u(p>) are plotted in fig. [2] for various width values. 

With u(<p) and f(z) defined by formulas (|2^|) and f )25|) . our patch model (122]) becomes: 



V(p> — p*,z — z*) = V(p — z — z*\ Az, Ap) = K v(p — tp*,z — z*] Az, Ap) 



(28) 



V* 



1 + cos (ip — p*) 



1 - 2A* cos (tp - p*) + A; 



exp 



V2Az~* 



The corresponding Fourier coefficients are found by the expressions ( 123]) . ( 124]) . and ( 1271) as 

2" 



V n (k) = 
Vo(k) 



1-A 2 

2vrA^V;^-^Al n| exp 



4A* 
1-A, 



exp 



/ kAz* 
(kAzX 



-i(nip t +kz t ) 
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±1,±2 



-ikz, 



(29) 
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Certainly, these functions are smooth enough to satisfy conditions (jl]), in fact, the function 
([28]) and all its derivatives in and z are squarely integrable, or else V belongs to the 
Sobolev space H p for any p > 0. The picture of equipotentials of the patch ( l2"gj) normalized 
by the maximum voltage is in fig. El 



V. SINGLE PATCH AT EACH OF THE ELECTRODES: A PICTURE OF PATCH 
INTERACTION 

We consider now two patches of the form (1281) . one at the inner, the other at the outer 
boundary. Patch voltages in the boundary conditions ([2]) become (i = 1, 2): 

G((p',z) = V(<f' - <fi,z -zx), H(<f,z) = V((p- <p 2 ,z-z 2 ) ; (30) 

_ ( z-Zj \ 2 
\V2AzJ \ ' 

where, according to the relation ( 1261) between A and the angular width, Aip, 
< Atpi < it, < Azi < oo, —n < Lpi < 7r, — oo < Zi < oo, % — 1, 2. The forces cor- 
responding to these distributions are derived in the appendix. We study here the interaction 
of two identical patches, Az x = Az 2 = Az, A</?i = A(p 2 = A<p, and V\ = ±V 2 = V . 



V(ip -(fi,z- z^ = Vi 



(1-Xif 1 + cos (ip - ifj) 
2 1 — 2Aj cos (ip — <fi) + A; 



exp 



A. Transverse Force 



1. Transverse force due to patch and uniform potential interaction 



By formula (1A1I) of the appendix, the transverse force due to the interaction between the 
uniform potential and patches reduces to the following expressions: 
Fi nt [WAz 



x 
1 



[cos 2tpi =f cos 2<y2 2 ) cos A<p 



[W Az . o a I / \ 
TT\ sin Aip< (C0S(y9i =FCOS(^2j — 

V 2 a [ 
2(1 T l) 



1 + cos Aip 



y 

— [(sin 2ipi =f sin 2<p 2 ) cos Aip] 



jpint 



[WAz . 2 f . . , 

-7tw — — sm A(/)< (siii(/3i =)= sin(^ 2 J 

i Zi (X I 



(31) 



— [(sin 2y?i =F sin 2<^ 2 ) cos A^?] — — 

tl tl 



2(1 Tl) 



1 + cos Aip 

Here F is the characteristic force defined as 



— (cos 2<pi =F cos 2ip 2 ) cos Aip 



F = e V V-(a/dY 



(32) 
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and V~ is the uniform voltage difference ([3]). The minus or plus sign is taken for the patch 
voltages of the same or opposite sign, respectively; the signs of the charges induced by patches 
on the cylinders are opposite in the first case, and same in the second. As expected, the 
magnitude of the transverse force is inversely proportional to (a/d) 2 . It is also proportional 
to the relative axial width, Az/a, and entirely independent of the axial positions of the 
patches. The dependence on the angular width is more complicated: the maximum force is 
at Aip = 7r/2, as prompted by geometry; for a small width the force goes to zero as (Aip) 2 . 
When Aip = it, i.e., the patches are the 'belts' of voltage uniform in ip, the zeroth order 
force vanishes, and the total becomes proportional to the shift and directed along it: 



This is zero when the patch voltages are equal: the forces from each of them have the same 
magnitude and opposite directions (see more on this below). 

The main contribution, i.e., the force in the centered position is best characterized by its 
polar components, namely: 



The total force is a superposition of the two forces from each of the patches. They act along 
the radial direction to the corresponding patch center, and the total is a vector sum of these 
two radial vectors, see figs. HI That is why the interaction force vanishes when the patches 
are one opposite the other (<pi = ip 2 ): their contributions, aligned and of the opposite signs, 
exactly cancel each other. So here a patch behaves as an effective point charge, q e ff, m the 
uniform radial field E u = —V~/d: the force due to it is just q e ffE u . The effective charges 
are readily found by comparison with the above expressions of the force. 

The forces of the first order are directionally coupled to the shifts, meaning an x-force 
depends on the y-shift, and vice versa. The forces consist of a constant term and the second 
harmonics of the patch angular position. 




(33) 
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2. Transverse force due to patch interaction 



Its general expressions (1A8|) and (lAlOj) simplify for the same size patches: 



M =F2>liexp 



Zi ~ Z 2 
2A^ 



(COS </?l + COS tf 2 ) 



X 

1 



2Ao + A/" 2 (cos 2<^>i + cos 2ip 2 ) =F 2 A/f + M 2 cos (</?i + ip 2 ) exp 



Zl ~ Z 2 

2Az 



y_ 
d 



Af 2 (sin 2ifi + sin 2<^ 2 ) =F 2A4 2 sin (<£>! + <£> 2 ) exp 
o 



z 
Az~ 



. \ z i — z 2 

M. i (cos (fi + cos ip 2) — -r — exp 



z\ - z 2 
2Az 



z l — z 2 

2Az 

2 



f- = 2(-) 3/2 



¥ sin2 (t) 



M =F 2A4iexp 



gi - Z 2 
2A;z 



(sin^i + siny2 2 ) - (34) 



x 

1 



J\f 2 (sin 2(pi + sin 2tp 2 ) =F 2A4 2 sin (<pi + <^ 2 ) exp 



zi - z 2 
2Az 



d 



2Ao — A/2 (cos 2^! + cos 2(p 2 ) =p 2 yM — M 2 cos + ip 2 ) j exp 





z 
Az~ 



K A f ■ ■ \ Z l ~ Z 2 

Mi (sin ifi + sin (p 2 ) — exp 



zi - z 2 
2A^ 



Zl - z 2 

2A2 

2 



Here the characteristic transverse force, F tr , is almost as in the previous case [formula (1321)]. 



(35) 



with just a natural replacement of V with Vq. The coefficients involved are found by 
combining the formulas (1A5|) . flA6j) and (1A7I) of the appendix with the representation (IA12I) : 
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(1 + 4A-3A 2 ) ; 



M 
A4i 
M 2 



1 - A 
4 

1-A 2 



1 + (1 + A) 



1-A(1 + A) 



2 cos(y?i - yj 2 ) - A 2 



2D 

1 + A 2 - 2cos(y?i - (p 2 ) 
2D 



D = 1 - 2A 2 cos (y?i - (p 2 ) + A 4 ; (36) 



1-A 2 f 1 + A 



A 



2cos(v?i -<p 2 ) + A(l + A) 



cos 2(<pi - ip 2 ) - A 2 cos(y?i - ip 2 ) 



D 
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As before, the signs =F correspond to the case of the same or opposite signs of the patch 
voltages. The force (|34"|) is again proportional to {a/d) 2 and Az/a. The dependence on the 
angular width here is even more complicated than in the previous case; still, the force is 
oc (Ay?) 2 in the narrow patch limit. For A(p = tt, an analog of the formula holds: 



F* = 2(nf 2 F tr (^) ^exp 



Zl - z 2 
2Az 



Po 
d 



This force, however, does not vanish for identical voltages, unless the patches are one right 
against the other. 

The nature of the force due to the patch interaction is most clearly seen in the main term 
corresponding to coaxial cylinders with no axial shift: 
a) V 1 = V 2 = Vp 

21 1 



5 



cos • 



(px - y? 2 



.Ml - 2.M 1 exp 



b) V 1 = -V 2 = V 



£=L = 4(7 r) 3 / 2 



Az . 2 {A<p\ 

— sin 

V 2 J 



cos ■ 



ifX - ip 2 



Af 1 + 2M 1 exp 



Zl - z 2 
2Az 



zi - z 2 



2Az 



The direction in the cylinder cross-section is along the bisectrix of the angle subtended 
by the patches, no matter what the voltage signs are:tan# p = tan[0.5(</?i + y^)]- Thus for 
ifi = tp 2 the total force is maximum because the contribution of each patch doubles, while 
for (pi — (p 2 + 7T the total force is apparently zero. The force is a combination of two terms: 
one is proportional to A/i depending only on the angular patch width and non-negative, by 
the formulas f )36|) . The other is the exponent of (zi — z 2 ) 2 with the coefficient Aii depending 
on \ip x — V^l- The coefficient M. Xl along with A^o; an d M 2 , are shown in fig. Inversus the 
angular distance. They have a sharp maximum at (fx — <p 2 , and drop quickly away from it; 
the maximum is the sharper, and the drop the faster, the smaller the width of the patch is. 
This is a strong manifestation of screening of patch charges: the patches interact strongest 
of all when their centers are right opposite each other; the interaction drops when one charge 
stops 'seeing' the other due to the obstruction by the inner cylinder. 

In fig. [Hl-F^ is plotted as a function of the angular distance \<px — tp 2 \ f° r both cases, a) 
and b). The two curves differ significantly for the moderate angular separations, and tend 
to zero when ipx — ¥2 ~^ tt- The dependence of F± on the axial distance \zx — z 2 \ is shown 
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in fig. [7J The force does not vanish at infinity, but goes instead to the asymptotic value 
common for the two cases. 

As prompted by the similarity between the coefficients Ai — Ai 2 and Af — Af 2 , being just 
some bounded functions of Atp, the terms in the force ( EM]) proportional to the transverse 
shifts have the structure, and thus the behavior, similar to that of the zeroth order expres- 
sions. The term with the axial shift is rather different, first of all because of an additional 
factor, (d/Az). Thus z° no longer compares to the gap, d: quite naturally, it is the ratio 
(z°/Az) that stands as a small parameter at this part of the force. The other peculiarity is 
the additional factor proportional to z\ — z 2 . This part of the force vanishes in both limits, 
[z\ — z 2 )/ Az — > oo and (z\ — z 2 )/ Az = 0, with the maximum magnitude at \zi — z 2 \ = y2 Az 
as seen in figJHl Finally, it is proportional to the first harmonics of the angular patch posi- 
tions and the coefficient M.\. The coefficient, and thus the force, is maximum when cp 1 = <£> 2 , 
then it decreases monotonically as \ipi — ip 2 \ increases, fig. [51 



B. Axial Force 



As demonstrated in section IHH the axial PE force is only due to the patch interaction. 
Its general expression ( 1A11I) reduces in our case to: 

21 ( 



^ = ±2^ sin 2 f^*-^ 



F f 

x 



ax 




V 2 J Az 

y° 



exp 



z\ - z 2 
2Az 



^ Mi (cosy?i + cosip 2 ) - -^-A^i (sin<£>i + sin<£> 2 ) 



Mo 

z° 
Zi - z 2 



(37) 



\-( 



z\ - z 2 



V V2A; 



The coefficients Aio, Ai\ are found in the formulas (1361) . and the characteristic force is 
inversely proportional to the relative gap, d/a, and not to its square, as before: 



Fax = e V a/d 



(38) 



The sign of the force (1371) switches from plus to minus for patch voltages of the same 
or opposite signs, respectively. The most striking feature of the axial force is its overall 
dependence on the ratio [z\ — z 2 ) / Az. As one expects intuitively, the axial force tends to 
zero in both limits Az — > and Az — > oo, due to the Gaussian exponent of the above 
argument multiplied by this same argument. It has the maximum 7r 3 / 2 e -2 sin 2 (Acp/2) F ax 
at \z\ — z 2 \ = \^2Az, same as the axial shift part of the transverse force from fig. [HI 
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The overall dependence of the axial force on Atp is just like that of the transverse force; 
particularly, for belt-like patches, Aip = tt, the force is: 

( Z\ — Z 2 \ * I , Z" _ ( Z\ — Z 2 ) 

exp 



EL = ± ^3/2 gi z Z 2 

F ax Az 



2Az 



1 



z\ - z 2 



1 



( Z1 ~ z 



Here the main term vanishes when the two patches are in the same cross-section, z\ — z 2 , but 
the 'correction' proportional to z° does not: F^/F ax = =f tt 3 / 2 (z°/Az), Aip = tt, z\ — z 2 . 

The first order force proportional to the axial shift has a factor {z\ — z 2 ) 2 , instead of 
{z\ — z 2 ), in front of the Gaussian. This contribution changes sign at z\ — z 2 = y2Az, and 
is maximum when the patches are in the same z plane, and it as illustrated by fig. [9j 

The axial force decreases with the angular distance between the patches. The zeroth order 
component of the axial force together with the first order term proportional to the axial shift 
depend on \<pi — (p 2 \ through the coefficient A4q only. The two other components due to 
the transverse shifts are proportional to the harmonics of the patches angular positions and 
the coefficient M.\\ they both vanish in just one case, when ipi — ip 2 = tt. Despite these 
minor differences, each of these terms follows basically the characteristic behavior of the 
coefficients A^o or M.\ versus the distance \y>\ — ip 2 \, as presented in fig. [51 



VI. ESTIMATES OF AXIAL PATCH EFFECT FORCE FOR STEP. 
CONCLUDING REMARKS 

A. Basics of the STEP Experimental Set-up and Some Requirements 

STEP, a medium size scientific satellite (< 900 Kg), will be put into a drag-free earth 
orbit at the altitude of ~ 550 Km. It is to measure the relative free fall acceleration of pairs 
of test masses (TM) of different materials with an accuracy of 10 _17 m/sec 2 , to determine the 
equivalence of the inertial and gravitational mass with an uncertainty 6 orders of magnitude 
smaller than the existing results 15], or find violations of the Equivalence Principle (EP; in 
other words, Universal Free Fall) somewhere between 1 part in 10 12 and 1 part in 10 18 . 

STEP will fly four differential accelerometers (DAC), each with a pair of TMs shaped as 
coaxial cylindrical shells. The cross-section of the DAC is shown in fig. [10J An electromag- 
netic system of magnetic bearings and capacitors keeps TMs alligned and centred to within 
< lnm. The transverse degrees of freedom are constrained, the two axial ones are left free, 
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so the axial motion and the rotation about the TM axis are most important. So we discuss 
only the axial PE force; the axial torque will be examined in the final part of this paper. 

The readout of the differential longitudinal displacement for each TM pair is provided by 
a Superconducting Quantum Interference Device (SQUID), sensitive to an acceleration of 
3 x 10~ 18 m/sec 2 in an accumulation time of at least 20 orbits, i.e., ~ 2 days. Nominally the 
DACs are kept inertially fixed and the EP violation signal will be at the orbital frequency, 
f or b = 1.74 x 10~ 4 Hz. However, changing the signal frequency from one measurement 
session to the other helps to discover and remove systematic readout errors. So modulation 
of the frequency is planned by rolling the satellite (and the DACs) about the normal to the 
orbital plane at some frequency f ro u > 2f ar b (the DAC axes lie in the orbital plane). With 
this procedure, the science signal will be at the frequency f s = f ro u ± / or & bounded from 
below as 

f a > 1.74 x W~ A Hz . (39) 

The STEP design rejects all perturbations able to mimic or mask the signal from the 
TM free fall at the level of the target accuracy. Thus the magnitude of any perturbing axial 
acceleration at the signal frequency should be at least less than 10~ 17 m/sec 2 . In addition, 
axial motion of TMs is allowed only within certain range, as required by drag-free control 
and SQUID readout limits. The total axial shift must not exceed some limit, z^ ax , whose 
value may be adjusted while developing the instrument; currently it is z® max = 1 \im. 

The SQUID sensor will keep the TM position in the body by applying a restoring (spring) 
force to counteract any D.C axial force, such as the one due to patches. Its control authority 
in terms of the range of the corresponding oscillation period is 300 — 1000 s. For a TM of 1 Kg, 
this converts into the spring constant, k, within the range 3.9 x 10~ 5 — 4.4 x 10 -4 N/m. Its 
maximum value k max = 4.4 x 10 -4 N/m provides the maximum force that can be controlled 
according to F max = k max z^ nax = 4.4 x 10~ 10 N, or maximum acceleration 

dmax = F max /l Kg = 4.4 x 10- 10 m/s 2 . (40) 

The related criterion is: any distribution of patches is acceptable as soon as the axial D.C. 
acceleration due to them does not exceed a max . 

We consider a TM and its magnetic bearing as a reference case of our pair of cylinders 
(the gap between them is at least 3 times smaller - and the force equally larger - than the 
gap between the TM and electrodes, see fig. [10]). First we assume each cylinder to carry 
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just one patch of the same sizes and magnitudes; then we consider an arbitrary number of 
such patch pairs. Using slightly different notations the formula ( 157)) for the axial force with 



x° = y° = becomes: 



±7T 3/2 F ax sin 4 ]// cxp 



z\ - z 2 



H = /i(A, \(f! - (p 2 \) = 1 



(1 + A) 



2Az 

2 



Z\ ~ Z 2 



Az 



Az 

cos(ip 1 - ip 2 ) - A 2 



1 - 



z\ - z 2 



( V2Az 
A = cos Aip 



(41) 



D is 



2 1 - 2A 2 cos (<f i - (p 2 ) + A 4 ' 

the characteristic force F ax is defined by the expression (1551) . To 1. o. the force 
constant (unless the patches move on the surfaces, which has not been observed so far, or 
the cylinders rotate); the correction proportional to the axial shift, z°, can produce harmonic 
oscillations near a stable equilibrium, or exponential runaway from an unstable one. We take 
the patch voltage Vq = 10 mV for our estimates. This seems a plausible number for patches 
at low temperatures (see some relevant result for the GP-B experiment in [TJ). The TM 
parameters used are: TM radius a = 2.3 cm (outer TM; the inner radius, and hence the 
force, is about 5 times smaller), TM height 2L = 14 cm, the TM to magnetic bearing gap 
d = 0.3 mm (d/a fa 10~ 2 ). The TM masses vary from 0.3 K g to 2AKg, so to make the 
rescaling easy, we give all the accelerations (or specific forces) per 1 Kg of mass. 



B. Constant Axial Acceleration due to Patches 

By the expression (l4Tj) . the constant part of the axial patch effect force is 

F z = ±27r 3 / 2 F a ,sin 4 (^f)t*z exp (-z 2 ) , z = (z, - z 2 )/2Az ; (42) 

the ballpark number for the resulting acceleration comes from the characteristic force (155)) : 

a ax = F ax /1 Kg = e V 2 (a/d)/l Kg w 6.8 x 10" 14 m/s 2 . (43) 

It is almost 4 orders of magnitude smaller than the one required by the formula ( 14*01) . Next is 
an upper bound on the acceleration from the force (142D with two last factors set at their max- 
imum. Since the maximum of \i with regards to | f\ — f 2 \ is 1 + sin 2 (Aip/2) /2 sin 2 (Aip/2), 
and the last factor's maximum is l/v2e when z = l/y/2, the bound is 



a z < 



'27T 3 



g a ax snr y—j ~ 3 -2 x 10 



-13 • 2 f A( f\ I 2 
sin I "g - ) m / s 



(44) 
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It is still more than 3 orders of magnitude smaller than the required value ( HO"]) for any Aip, 
and drops as its square when it becomes small. 

We now consider the mean value of the force (I42p over all angular \ipi — ip 2 \ and the axial, 
z patch distances (Zo = L/2Az). The average in the angle is unity while that in the distance 
is (/ )- 1 /o°5exp(-5 2 )^ = (Zo)- 1 [1 -expHo) 2 ], so 

a z = 2^ 2 a ax ^- [l - e"^ 2 ^ 2 ] sin 4 < 7.6 x 1(T 13 ^ sin 4 (^f) m/s 2 . (45) 

Two new pleasant features here as compared to the estimate (jHJ) are: the fourth power 
instead of the square of sine, and a factor (Az/L), which is less than unity in any case, and 
is expected to be essentially smaller. The upper bound (j4"5|) satisfies the condition PUj) with 
the margin of 3 orders of magnitude for any patch sizes. 

However, this is hardly the case with some number, N, of patches pairs identical in sizes 
and magnitude, since N may be quite large. Those pairs generate N 2 interactions each 
giving the force ( 142]) . with either sign. Assuming the signs of the patches random, single 
contributions do not cancel each other completely, but sum up to a factor \A/V 2 = N for 
N ^> 1. Thus the total force is reasonably estimated as iV times the average of the force 
( H2]) . The resulting acceleration satisfies the condition ( 1401 if the following inequality is true: 

jV^siW^ <5.8xl0 2 . (46) 

For N up to 600 the condition is satisfied by the patches of any sizes. Limitations on the 
latter start with N ~ 1000, they are not too restrictive until N ~ 10, 000: even for this 
number and the angular size A(p = 60°, any axial size of the patch is still acceptable. 

Moreover, the upper bound for N is available through the patch sizes, if the patches do not 
overlap within their nominal widths 2Az, 2Atp. The effective patch area is (2Az) x (2aA<p), 
the total surface area of the cylinder is (27ra)(2L), so iV is bounded by the ratio 

, T „ AnaL 7r L 

N < = . (47) 

~ 4(aAip)Az Atp Az v ' 

Introducing this to the inequality (j4"6"]l . we obtain a universal estimate in terms of A(p only: 

J-si^(^)< 180, (48) 

which, of course, is always true. So, non-overlapping patches covering the surfaces completely 
satisfy the STEP limitation on the DC acceleration, if the averaged force is used and the 
assumption of random signs of the patch potentials is valid. 
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C. Harmonic Oscillations and Exponential Runaway due to Patches 



The zero order force (14T|) vanishes at z\ = z 2 . The first correction to it, SF Z , is: 

5F Z = T vr 3/2 ^ sin 4 (^)/, (l - 2~z 2 ) exp (-5 2 ) , z={z x - z 2 )/2Az . (49) 

As explained in sec. IV B\ its overall sign depends on the signs of the patches and the distance 
between them, because the factor in brackets is positive or negative depending on whether 
z is smaller or larger than l/y2. So we denote u 2 = \SF X \ /l Kg, then the equation of TM 
motion in the axial direction near the equilibrium becomes: 

z'o = Tcu 2 z o } (50) 

and describes small oscillations for the minus sign (restoring force) , and exponential runaway 
otherwise. In the first case, the oscillation frequency is / = oj/2-k > 0, in the second the 
TM drifts away exponentially with the characteristic time r = 1/oj. 

The maximum of u, needed for an upper bound for the frequency / or a lower bound for 
the time constant r, is clearly attained when z — (pi — (f2 — 0, so: 



u < sin 2 (^j = 6.2 x 10" 7 ^= sin (^j rad/s , (51) 

with Az in meters. Remarkably, this estimate of u decreases when Aip goes down, but 
increases with the decrease of Az. The frequency of TM oscillations should be below the 
minimum science signal frequency (l39j) . 

/ < 1.74 x 10~* Hz , (52) 

to avoid the signal corruption. Based on the estimate (jBTj) . this condition holds when 

m . (53) 



Az > 3.2 x l(r 7 sin 2 (~ipj 



Even if sin (A<p/2) = 1 (circular patches), the axial width should be only about two hundreds 
nanometers to satisfy this. 

In the runaway case the characteristic time must satisfy 

r»T o6s = 1.7xl0 5 s, (54) 
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to avoid saturation: if it is true, then any initial shift smaller than z max = 1 jjm will not 
grow critically during a science session. This criterion and the estimate ( 151]) lead to 



Az > 1.1 x l(T 2 sin 2 m 



(55) 



five orders of magnitude more restrictive than the bound (|53|) . However, a small enough 
angular size helps to meet the requirement fl54|) . The maximum acceleration allowed in the 
runaway case when the condition (I541) holds is implied by the estimate (!5T|) with Az min 
taken from the r.h.s of the estimate (155]) and the maximum shift z^x = 1 l^rn: 

sin 2 (A<p/2) 



LO Z— 



3.8 x 10 
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AZr, 



m/s 2 = 3.4 x 10~ 17 m/s 



(56) 



Its value is 7 orders of magnitude smaller than the STEP requirement ( 140]) . 

Other estimates for the patch parameters are produced by the approach of the previous 
section using the expression of \5F Z \ averaged over the axial and the angular distances 
between the patches. As we know, the average of \x in the formula ( 149]) is unity, and 
the average over the axial distance is Jq° dz (1 — 2z 2 ) exp [— z 2 } = exp(— Z 2 ,) (as before, 
la = L/2Az). Thus the mean angular frequency is found to be 
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sin exp 
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(57) 



\2V2A Z/ 

with Az on the far right in meters, as usual. As compared to the estimate (|5T]) . here is 
an extra power of the sine of Atp, and, most important, the Gaussian exponent that drops 
sharply with L/Az growing. Due to this, the conditions ( 152]) and ( l54j) on the oscillation 
frequency and runaway time are satisfied for the patches of any sizes, by the estimate (157]) . 
Something close to this happens if we consider N pairs of patches instead of one. Let us 
examine the condition fl54"]) . which is more stringent than ( 152]) . We can again estimate the 
overall contribution of iV pairs of patches with random voltage signs as a/ N 2 5F z , and take 
the mimimum \/L(e/2) 1//4 for Az = L/\2. This gives 



N 1/2 sin 2 (Ay?/2) < 9.4 VZ(e/2) 1/4 = 2.7 ; 



(58) 



with L = 7 cm cited above. The inequality does not hold always, somewhat limiting the 
angular patch size to be small enough given N. However, assuming the patches do not 
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overlap, we can use the upper bound (|47|) of N with the above value Az = Lj ' \f2 in it. This 
turns the inequality (JHE} into sin 2 (A<^/2)(A^/2) _1/2 < 1.8. The maximum of the l.h.s here 
is always less than unity, so the condition (158]) holds always for non-overlapping patches. 

Concluding this section we remind that all the numbers in the estimates were obtained 
using the patch voltage Vq = 10 mV, the TM mass 1 Kg, and the TM radius-to-gap ra- 
tio a/d = 77. All the accelerations scale as voltage square, proportional to this ratio, and 
inversely proportional to the mass; the frequencies and inverse runaway times scale propor- 
tional to the voltage, Jajd, and the inverse of the square root of the mass. 

D. Concluding Remarks. Perspectives of Patch Force Modeling 

The main message of the above estimates is that the STEP requirements f)4Up . (I52p and 
(I54j) can be met by patches whose sizes and number are appropriately limited; the limits 
seem not too restrictive and practically achievable. However, those estimates were derived 
under a number of simplifying assumptions. Perhaps the strongest of them was that all the 
patches differ at most by the sign of the voltage, having the same voltage magnitude and 
spot sizes. The results obtained in this paper allow, in fact, for an essentially more realistic 
patch force modeling, which can be done in the following way. 

One represents both patch potentials V a ((p ,z) and V&(<£>, z) as a superposition of some 
number, AT 0) &, of model patches ff28l . 

V^tp, *) = £ V(<p ~rt,z- z£; Az%, A^ n ) = 

n=0 

E^«(WS,2-* A<,A^), i* = a,b, (59) 

n=0 

with different voltages, sizes, and positions. By the formulas of section II III the patch effect 
forces corresponding to the distributions (J59~j) can be explicitly calculated, as it is done here 
for a single patch at each of the boundaries. Being cumbersome, this general calculation 
is otherwise straightforward, without any new technical difficulties. It leads to the force 
expression as a quadratic form of the patch voltages V£, with the coefficients depending on 
all other parameters in a known way. 

Having these formulas at hand, one then carries out simulations by specifying parameter 
sets in various ways and computing the patch forces. One can pick the parameters randomly, 
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and eventually come up with the patch force statistics. One can also use any lab information 
on the patch distributions, arranging for a semi-random patch sets, as was done, for instance, 
when simulating magnetic trapped flux distribution on GP-B rotors 16|. Such exhaustive 
analysis can be strongly recommended before the STEP flight. On the other hand, the same 
general formulas for the transverse forces can be used for fitting control effort data obtained 
during the experiment, for restoring the voltage patch patterns on the proof masses and 
bearings. Once the latter are known, the axial forces can be computed, and the systematic 
experimental error due to them can thus be bounded. Of course, all this is applicable to 
any experimental set-up with cylindrical geometry. 
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Appendix A: Calculation of the Transverse and Axial Force for a Single Patch at 
Each of the Cylinders 



Here we consider one patch at each boundary described by our patch model ([28]) . i.e., for 
the voltage distributions (130|) . The force ( |T8|) due to the interaction between the uniform 
potential and patches is a linear function of the Fourier coefficients of the boundary patch 
voltages computed for k = and n = 0,1,2. Thus formulas (12"§|) allow one to obtain 
(h = s/2A Zi , % = 1, 2): 
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2 x / »» * \ '^/i 2 

The transverse force due to the patch interaction is essentially more cumbersome to derive, 
with the additional difficulty of computing some integrals in k and sums over n containing 
products of Fourier coefficients. We start with F£ from the first of the formulas ( !2~TI) . Using 
the expressions ([231) and (13"Uj) . combined with the equality (1211) . we find: 
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Here we have introduced the following notations for the coefficients: 

1 00 

M q = M q (ipi, Xi,(f2, A 2 ) = — £ M „(A 1 )e-^ 1 < +(? (A 2 )e 



i(n+g)ip2 . 



(A3) 



N q (X) 



2tt 



^ u n (A)< +? (A) = M(0,A;0,A); g = 0,l,2 



(A4) 



u n (A) defined by the formulas (|2"7|) . Coefficients (1A3I) are clearly symmetric in the pairs of 
arguments (p, X corresponding to each of the patches, M q ((pi, Ai, <p 2 , A 2 ) = M q (ip 2 , A 2 , <pi, Ai), 
which we have already used in the above expressions. Now, formulas ( I27I) lead to the explicit 
sums of the series (1A3|) . since they reduce to geometric progressions: 
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M 2 ((fx, Xi,(f2, A 2 ) = 
e 2 ^ 2 (l + A 2 )A 2 



(1 - AQ (1 - A 2 ) 
8 



2 v ; D 
(1 + Ai)(l + A 2 ) c ^ 2) l 



D = 1 - 2AiA 2 cos(^i - <p 2 ) + (AiA 2 ) 2 . (A6) 
By the definition (1A4[) . coefficients N q (\) are found as a particular case of the expressions 



1A5I) when Ai = A 2 = A and c^i = c^ 2 = 0: 

A (A) = ^(3-A); iV^A) = i-^(2-A) ; A 2 (A) = i-^!(l+4A-3A 2 ) . (A7) 
8 8 Id 



For the closed form of Fg, two well known integrals in k are used: 
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In the case z\ = z 2 and l\ = l 2 = I the first integral becomes: dk exp 
Thanks to all the above results for the series and integrals, we obtain the expression for F£, 
to 1. o. in the parameter r°/d <C 1: 
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Calculation of F£ is pretty similar, and its result, to linear order in r°/d <C 1, is: 
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Finally, we find a closed-form expression for the axial force. Starting from the last of 
the formulas ( 12~Tj) . we just need to repeat the same steps as with Fg. The only significant 
difference is that here we need a slightly different integral, 
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In this way, using definitions ( 1A9I) and ( 1A5I) . we arrive at the final expression for F£: 
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The calculation of the patch effect forces for the case when a single patch is placed on 
each of the cylinders is now completed. The expressions (1A8|) . (lAlOj) and (1A11I) are for the 
patches with different sizes and magnitudes. When those are identical, the results simplify 
essentially due to Ai = A 2 = A, l\ = l 2 = I, \Vi\ = |V 2 |, in particular: 

iVo(A) = (1 - A) Af ; A^A) = (1 - A)M ; iV 2 (A) = (1 - A) A/ 2 ; 

M (A,^i,^ 2 ) = (l-A)wM ; (A12) 
Mi (A, <f i, if 2 ) = (1 - A) .Mi (A, fx - if 2 ) [(cosy?i + cosf 2 ) + i (sin</?i + sin^ 2 )] ; 
M 2 (A, </?i, ^2) = (1 - \)M 2 (\ifi - f 2 ) [cos (fx + f 2 ) + «sin (if i + if 2 )\ . 



The coefficients Mi and Aii here are given explicitly by the formulas (j36p . 
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FIG. 1: Geometry of the problem and coordinate systems 
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FIG. 2: Azimuthal profile of the patch for various Aip = 1°, 70°, 90°, 110°, 179° 
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FIG. 3: Equipotentials of the patch model for Aip = n/8 




FIG. 4: Interaction force due to two patch with a) same sign of the voltages; b) opposite sign of 
the voltages 
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FIG. 7: Normalized transverse force vs. the axial separation of the patches 
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FIG. 9: Normalized z-component of the force due to the axial shift as a function of the axial 
distance between the patches 




a = 5xlO- 3 



FIG. 10: Cross-section of the STEP differential accelerometer 
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